Longitudinal Broadening of Quenched Jets in Turbulent Color Fields 
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The near-side distribution of particles at intermediate transverse momentum, associated with 
a high momentum trigger hadron produced in a high energy heavy-ion collision, is broadened in 
rapidity compared with the jet cone. This broadened distribution is thought to contain the energy 
lost by the progenitor parton of the trigger hadron. We show that the broadening can be explained 
as the final-state deflection of the gluons radiated from the hard parton inside the medium by soft, 
transversely oriented, turbulent color fields that arise in the presence of plasma instabilities. The 
magnitude of the effect is found to grow with medium size and density and diminish with increasing 
energy of the associated hadron. 



The emission of hadrons with large transverse momen- 
tum is observed to be strongly suppressed in central col- 
lisions of heavy nuclei at the Relativistic Heavy Ion Col- 
lider (RHIC) [2,i]. The suppression is understood to be 
caused by final-state rescattering of the leading parton 
in the dense medium produced in such collisions, causing 
it to lose energy. Within the framework of perturbative 
QCD, the leading process of energy loss of a fast parton 
is gluon radiation induced by clastic collisions with color 
charges in the quasi-thermal medium 0, 0, Sj- Elastic 
collisions not followed by radiation may also contribute 
to the energy loss [fj. 

In this manuscript, we study the fate of the energy 
lost by the leading parton of the trigger jet within the 
medium. Measurements by the STAR collaboration have 
shown that the leading hadron is accompanied by a cone- 
shaped pattern of secondary hadrons, which closely re- 
sembles the jet cone around a hard scattered parton in 
nucleon-nucleon collisions @, Hj]. We will refer to this 
as the contribution from vacuum fragmentation (or "jet" 
yield). In addition, the trigger parton is observed to be 
accompanied on the near side by a wide ridge (or 
pedestal [lfj) of secondary hadrons, which extends along 
the beam direction over more than one unit of rapidity in 
both directions. We refer to this as the contribution from 
in- medium fragmentation (or "ridge" yield). The energy 
flow in the ridge-shaped structure grows with participant 
number and is concentrated in hadrons with transverse 
momenta pr < 2 GeV/c. 

The inclusive yield of associated hadrons receives con- 
tributions from the fragmentation of the leading parton 
as well as from the radiated gluons . The yield of asso- 
ciated hadrons produced in the fragmentation of the hard 
parton after its exit from the medium will be assumed to 
form the jet yield [12 ]. This naturally explains why the 
azimuthal and rapidity distribution of this contribution 
is identical to that from a jet of the same final energy 
produced in a p + p collision. Hadrons produced in the 
fragmentation of the gluons radiated within the medium 
or after absorption of those gluons will be assumed to 
form the ridge yield. In either case, the ridge shape will 
reflect the kinematic distribution of the radiated gluons 



at the end of their passage through the medium. 

It has been argued that the properties of the ridge 
must be associated with the longitudinal expansion of 
the medium (ToL fl3j ]. However, the precise mechanism 
by which the ridge acquires its peculiar shape has not 
been quantitatively explained. In a recent article, Ro- 
matschke [14| has shown that the momentum of a trans- 
versely propagating heavy quark is preferentially broad- 
ened in the longitudinal direction by elastic collisions 
in an expanding medium, if the expansion leads to a 
large anisotropy of the momentum distribution of the 
partons composing the medium. He suggested that this 
phenomenon could explain the ridge. We note, how- 
ever, that elastic scattering is probably not the primary 
source of energy loss for the light partons which initi- 
ate the observed ridge. We also note that the momen- 
tum anisotropy parameter £ of the expanding medium 



is related to its shear viscosity rj by the relation [15] 
£ ps 10 t]/(stT), where s denotes the entropy density, 
T the temperature, and r the time after the initial hard 
collision. For times relevant to the propagation of the 
jet through the medium (t > 1 fm/c), the large values 
of £ ~ 10 required to explain the extent of the ridge are 
difficult to reconcile with the low shear viscosity of the 
medium (rj/s < 0.3) required by the data 14. 1 X 61] - 



Here, we propose a different mechanism for the ori- 
gin of the observed ridge, also related to the longitudinal 
expansion of the medium. It builds on the recent in- 
sight that extended color fields are dynamically formed 
in the expanding medium due to the presence of plasma 
instabilities. Such instabilities have been shown to exist 
in any quark-gluon plasma with an anisotropic momen- 



tum distribution [18|, [19|, |20| . The nearly boost invari- 



ant longitudinal expansion of the matter formed in rela- 
tivistic heavy ion collisions necessarily induces an oblate 
momentum distribution of partons with respect to the 
beam axis. The plasma instabilities lead to the exponen- 
tial growth of soft modes of the glue field, which ulti- 
mately saturate due to the nonlinear self-interactions of 
the Yang-Mills field, resulting in a turbulent state of the 
quasi-thermal quark-gluon plasma [2l|, [22] • Such turbu- 
lent color fields currently present a successful mechanism 
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for the fast equilibration of the medium 18, 20] and serve 
as a source of an anomalous shear viscosity [15j ] . This nat- 
urally explains the ideal fluid behavior of the medium. 
Here we explore the effect of turbulent color fields on the 
gluons radiated by a hard parton in the medium. 

When the hard parton scatters transverse (in the 
medium co-moving frame) to the beam direction, its ac- 
companying halo of soft radiated gluons are deflected 
by the turbulent color fields as the parton traverses the 
medium. The deflection pattern is not isotropic around 
the jet axis, because the direction of polarization of the 
instability driven color fields is not totally random. De- 
tailed studies of the instability pattern [19( have shown 
that color magnetic field modes polarized transversely 
to the beam direction exhibit the largest growth rates 
and, therefore, dominate the field configurations at early 
times. The dominance of transversely polarized color- 
magnetic fields results in the preferential deflection of 
transversely propagating partons in the direction of the 
beam axis, with a deflection angle that changes inversely 
with the parton momentum. The resulting pattern will 
be one in which gluons of moderate pt from the jet cone 
will fan out along the beam axis to form a ridge as it is 
observed in the RHIC experiments. 

In order to be deflected by the soft color fields, the 
gluon must be radiated within the medium. Denoting 
the gluon's momentum along the jet axis by px and its 
momentum components perpendicular to the jet by 1^, 
the formation time of the gluon is t{ ~ pr/l±- Therefore, 
the fields will most strongly affect those gluons which 
carry a small fraction of the jet's longitudinal momen- 
tum (with short formation times) and are emitted at a 
substantial angle with respect to the jet axis. These are 
precisely the characteristics of the gluons radiated due to 
the scattering of the original parton inside the medium. 

Consider a hard scattering event with induced radia- 
tion as sketched in Fig. [1] A jet initiating parton is pro- 
duced when an incoming light-like parton is struck by 
a virtual state with forward energy E and virtuality Q 2 . 
The final state parton may be produced off-shell by up to 
Q 2 . In cases where it is off-shell, it may reduce its virtual- 
ity by radiating a gluon immediately. It may be produced 
closer to its on-shell condition and re-scatter multiple 
times prior to radiating the gluon. Following the radia- 
tion, both remnants of the original parton will multiply 
scatter off the soft color fields present in the medium. 
The amplitudes for multiple scattering of the hard par- 
ton before and after the emission of the gluon, as well as 
the scattering of the radiated gluon itself will lead to a 
destructive interference in the forward collinear region 
of the gluon phase space, the Landau-Pomeranchuck- 
Migdal (LPM) effect. This interference suppresses the 
production of very collinear radiation with long forma- 
tion lengths. As a result the final radiation is produced 
over a short distance and then decoheres from the ini- 
tiating parton. The yield of inclusive gluon radiation in 



this setting has been studied extensively 0, [B|, [H, 0, . 
Here, we are interested in the more exclusive probability 
of emission of a gluon with large transverse momentum 
relative to the parent parton. 

To estimate this distribution, we adopt a simplified 
factorized approach illustrated in Fig. [TJ The short dis- 
tance process which describes the production of the par- 
ton and its multiple scattering in the dense matter fol- 
lowed by emission of the hard gluon is represented by 
the dark blob. It includes the two processes commonly 
referred to as hard-hard and hard-soft double scatter- 
ing [2|| and whose interference leads to the LPM ef- 
fect for very collinear radiation 25 , 2(| . The subsequent 



soft scattering of the gluon and the parent parton in the 
medium is treated as independent from the production 
and radiation process. The factorization is indicated by 
the rectangular boxes in Fig. [TJ 

The remnant jet and the radiated hard gluon then at- 
tempt to exit the medium prior to fragmenting into a 
jet of hadrons. On their way out, these may still endure 
multiple soft scattering off the dense gluonic field of the 
medium. Such space-like exchanges are most often too 
soft to induce further radiation, however, the cumula- 
tive effect of many such scatterings may lead to a non- 
negligible deflection away from the direction of propaga- 
tion. While such further re-scattering has a diminished 
role in the modification of the fragmentation function of 
the parton 27], its additive effect on the relative trans- 
verse momentum of the radiated gluon with respect to 
the leading parton is more pronounced. 



FIG. 1: A hard scattering followed by multiple soft re- 
scattering diagram 

To obtain the triggered gluon distribution before the 
final-state broadening, we take the ratio of the cross sec- 
tion a qg for radiating a gluon with momentum (pt 2 ,1_l) 
in the short distance process depicted in the left-hand 
square in Fig. [TJ and the inclusive cross section a q for 
production of quark with transverse momentum p^ and 
rapidity y . This differential triggered distribution at a 
fixed impact parameter b and transverse location r is 
given as Ref . [28| , 



d 2 a qg ( PT2 ,l ± ) PT )^±±- 
da„ UlPTl ' PT2j 2Wl 



(1) 



x/ d(p(r + nO^[2-2cos(7 lL ()}. 
s> 
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In the above equation, £ is the distance traveled by the 
produced partons from the primary vertex along the jet 
axis h prior to the scattering in the medium. The max- 
imum allowed value Cmax for £ is the distance from r to 
the surface. The factor tjl = l\/ (2pT 2 Z\) is the inverse of 
the formation time of a gluon with momentum lj_ trans- 
verse to the jet axis. The final quark momentum fraction 
is z\ — Ptx I (jPTi + Pt 2 ) ■ The overall constant C (p^ , pr 2 ) 
accounts for the fact that the final parton momentum in 
the numerator and denominator of ([1]) is different. 

The gluon distribution of Eq. (Q} is now used as the 
starting point f(t — to) for the time evolution of the cone 
of radiated gluons around the leading parton due to the 
influence of turbulent color fields. Note, that this dis- 
tribution includes solely those gluons whose production 
is induced by scattering in the medium and does not in- 
clude any vacuum contribution to the radiation cone. In 
order to calculate the deflection of the gluons by turbu- 
lent color fields, we employ a Fokker-Plank equation: 



d_ 

dt 



v • V r 



V p D(p,*)V, 



f = C[f], 



(2) 



with the average parton phase space distribution 
/(p, r, t) and the diffusion tensor 15|, |29| 



Aj = 



dt'iFMnmm,*))- 



(3) 



Here F = gQ a (E a + v x B a ) is the color Lorentz force 
generated by the turbulent color fields, and C[f] denotes 
the collision term. is related to the usual transport 
coefficient for radiative energy loss (jet quenching coeffi- 
cient) q via: 



(4) 



The above equation may be interpreted as a generaliza- 
tion of the scalar transport coefficient q to a tensor q a p. 
The tensorial form does not influence the total energy 
lost. However, it does influence the angular distribution 
of emitted radiation. Such an effect is not included in 
Eq. (fTJ) and will be treated in a forthcoming publication. 
The current effort explores the effect of the anisotropy of 
Dij on the radiated gluon after its decoherence from the 
leading parton. 

For random transversely polarized color-magnetic 
fields, the diffusion term may be expressed in a simplified 
form El: 



V pJ D(p)V p 



- g 2 Q 2 (j? 2 )r m 
2(AT| - l)Ej 



-(p) 



l(p) 



(p) 



(5) 



where N c — 3 is the number of colors, the index _L 
denotes the vector components transverse to the beam 
axis, and r m is the autocorrelation length (or time) of 
the color-magnetic field along the trajectory of the par- 
ton. E p is the energy of a parton with momentum p and 



v = p/Ep its velocity. LW = ip x V p denotes the 
generator of rotations in momentum space. 

In the case of interest here, the distribution f(p, t ) 
represents the gluons in the radiation cone of the pri- 
mary hard parton, which have already decohered from 
the parent parton, as described above. Given the form of 
the diffusion term, the evolution equation is best solved 
by decomposing the momentum distribution function in 
terms of spherical harmonics: 



/(p,t)=X>m(p,t)y£(0p,^p). 



(6) 
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In the above equation, m denote the total and z- 
component of the angular momentum in momentum 
space. The evolution equations of the various compo- 
nents ai m , obtained by substituting ([6]) into ([2]), is given 
by 



a lm (p,t) = aUpM)e^ {p)[l{l+1) ~ m2]{t ~^. 



(7) 



Only the spherically symmetric mode I = m = is un- 
affected by the evolution, while all other modes decay 
in time. Among these, modes with m = for a given 
I decay faster than those with m ^ 0. As a result, the 
distribution begins to broaden in the z-direction (beam 
direction). The evolution time is limited by the depth 
at which the partons were produced and is thus finite. 
Partons produced at greater depth therefore experience 
greater broadening along the beam direction than those 
produced closer to the surface. 

The results of our calculations are presented in Fig. [2] 
These show the unnormalized associated gluon distribu- 
tion around a jet initiating quark as a function of rapidity 
T] for a fixed azimuthal angle = 0, and as a function of 
4> at fixed rj = 0. The dashed (dotted) lines show the in- 
put distribution /(to) obtained from expression (JTJ) ; the 
solid lines show the modified distributions at the mo- 
ment when the gluons leave the medium. Both plots 
correspond to the original jet vertex formed at a depth 
of 3 fm. The plots in the top part of Fig. [5] correspond 
to an original jet transverse momentum px t + pt 2 = 10 
GeV/c, while those in the lower part correspond to an 
original jet momentum of px x + Pt 2 — 20 GeV. In both 
cases, the gluon carries a forward momentum fraction of 
22 = 0.4. We assumed in our calculation that the soft 
gluon density in the medium varies with time in accor- 
dance with the boost invariant longitudinal expansion of 
an ideal ultrarelativistic liquid. The initial value of the 
transport coefficient associated with the intensity of tur- 
bulent color fields was taken to be q = 2.2 GeV 2 /fm. As 
is clearly visible in Fig. [21 the associated gluon distribu- 
tion is considerably broadened in rj but not in <p. The 
extent of the broadening depends on the energy of the 
parton and drops with increasing energy. 

In summary, we have shown that near side energy loss 
in terms of radiated gluons in the presence of turbu- 
lent color fields can explain the experimentally observed 
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longitudinal broadening of jet cones in heavy-ion colli- 
sions. Our model predicts that this broadening decreases 
with increasing energy of the associated parton (the ra- 
diated gluon which fragments to produce the associated 
hadron). An extraction of the transport coefficient q from 
a measurement of such broadening in experiment will de- 
pend on the energy of the radiated gluon and, by exten- 
sion, on the hadronization mechanism. The momentum 
range of the partons contained in the rapidity broadened 
jet cone is in the regime where hadronization occurs dom- 
inantly via parton recombination. We therefore expect 
an enhanced baryon to meson ratio of hadrons contained 
in the ridge. 



q(x=lfm) = 2.2GeV 2 /fm 




FIG. 2: The gluon distribution around the trigger quark in r\ 
at <j> — and in <f> at n = 0. The dotted and dashed lines show 
the initial distributions obtained after the gluons have been 
radiated due to interactions of the hard scattered parton with 
the medium. The solid lines show the distributions at the mo- 
ment when the gluons escape from the medium. The curves 
in the top frame correspond to an original jet momentum of 
10 GeV/c and a final gluon momentum fraction of zi = 0.4. 
Those in the lower frame correspond to a jet momentum of 
20 GeV/c and the same 22. 

As mentioned earlier, turbulent color fields generated 
by Weibel instabilities in the early quark-gluon plasma 
have recently been proposed as a mechanism for the early 
thermalization of the matter produced at RHIC and as 
the source of a small anomalous viscosity. The current 



work has identified the broadening of triggered jet cones 
in rapidity as a more directly observable manifestation 
of the presence of turbulent color fields in the matter. If 
our interpretation can be confirmed by additional mea- 
surements, it will demonstrate the existence of plasma 
instabilities generating these turbulent fields and estab- 
lish their role in the transport of hard probes through a 
quark-gluon plasma. 
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